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Q\ • Abstract 

m ! 

, A versatile method is described for the practical computation of the ex- 

' act discrete Fourier transforms (DFT), both the direct and the inverse ones, 

■ of a continuous function g given by its values gj at the points of a uniform 
Q , grid Fjy generated by conjugacy classes of elements of finite adjoint order N 

I in the fundamental region F of compact semisimple Lie groups. The present 

0^- implementation of the method is for the groups SU{2), when F is reduced to 

fH , a one-dimensional segment, and for SU{2) x SU{2) • • • x SU{2) in multidi- 

mensional cases. This simplest case turns out to be a version of the discrete 
cosine transform (DCT). Implementations, abbreviated as DGT for Discrete 
Group Transform, based on simple Lie groups of higher ranks, are considered 
separately. 

^ ' DCT is often considered to be simply a specific type of the standard DFT. 

■ Here we show that the DCT is very different from the standard DFT when the 
properties of the continuous extensions of the two inverse discrete transforms 
from the discrete tj € -F/v to all points t G F are studied. The following 
properties of the continuous extension of DCT (called CEDCT) are proven and 
exemplified. 

Like the standard DFT, the DCT also returns the exact values of {gj} on the 
A'^ -|- 1 points of the grid. However, unlike the continuous extension of the 
standard DFT, 

(a) the CEDCT function /Ar(t) closely approximates g{t) between the points of 
the grid as well. 

(b) For increasing N , the derivative of /Ar(t) converges to the derivative of g{t). 

(c) For CEDCT the principle of locality is valid. 

Finally we use the continuous extension of the 2-dimensional DCT, SU{2) x 
SU{2), to illustrate its potential for interpolation as well as for the data com- 
pression of 2D images. 
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1 Introduction 



The decomposition of functions integrable on a finite segment into Fourier series of 
trigonometric functions of one variable is a well known method (e.g. Pfl Ej) whose 
theoretical and practical aspects have been thoroughly investigated during the last 
two centuries in connection with its numerous applications in science and engineering. 
It is natural to question whether any attempt to add something to it is not in fact a 
reinvention of what has been found before. 

Our general goal, which goes beyond this paper, is to elaborate a new decompo- 
sition method of functions of n variables into Fourier series using orbit functions of 
compact semisimple Lie groups El , with the idea of (i) making it accessible to users 
who are not specialists in Lie theory, (ii) underlining the versatility of its practical 
implementations, and (iii) most importantly, demonstrating the fertility of the un- 
derlying theme of this approach. One can find complete answers to limited questions 
(like the values of a finite number of Fourier coefficients) replacing the Lie group by 
a suitably chosen set of its discrete elements. The choice of the discrete elements is 
clearly crucial. 

In the context of our goal, the SU (2) case results in the simplest example of the new 
method, even though that is the case where the potential advantages of the method 
based on the symmetry groups could be most limited. Remarkably, however, in 
this low dimensional space the discrete Fourier transform on the SU(2) group results 
in one type of famous discrete cosine transforms discovered in 1974 [5], or more 
exactly the DCT-1, according to the currently accepted classification (see [3 EI). 
Its comparison with the standard method is both misleading and revealing. It is 
misleading because it is considerably similar (see jTTj) to the standard Discrete Fourier 
transform abbreviated typically as DFT (e.g. [1113 GDI)- Nevertheless it is revealing, 
since it can help to understand better the underlying reasons why for many practical 
applications the DCT is proven significantly more efficient than the DFT. In this 
paper we consider the concept of continuous extension of the discrete transform, 
and show that the convergence properties of the continuous extension of the inverse 
DCT, abbreviated here simply as CEDCT, match very closely the properties of the 
canonical {continuous) Fourier transform (CFT) of smooth functions. Meanwhile, 
the continuous extension of the inverse DFT, abbreviated simply as CEDFT, does 
not result in a reasonable function at all. Note that for the sake of simplicity, if the 
"inverse" is not explicitly used, we will henceforth adopt DFT and DCT (or DGT in 
a more general sense) abbreviations for both direct and inverse Fourier transforms. 

In Section 2 we present the basics of the Fourier analysis on the SU{2) group, 
which also demonstrates the general formalism used for the Fourier transforms on Lie 
groups. We show that in practice the Fourier transform of a class function of SU{2) 
to the orbit functions of this group is reduced to the decomposition of a discrete 
function {gj. \ k = 0, 1, ... , N} defined on the A^-interval grid of variable t G [0, Tq] 
onto the series of (A^+ 1) cosine functions (including cosO = 1) of the harmonic order 
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n — j/'^ < N/2. The basis for the DCT series is thus composed of the first N half- 
harmonics of the cos(27rt/To) function, i.e. cos(27rnjt/To) = cos(7rjt/To), meaning 
that the harmonic order of these functions may be formally both integer for j even, 
and half-integer for odd j. This approach is compared to the standard method of DFT 
where the given {gk} is decomposed into the trigonometric polynomials of 6^*^'^"*/-^°^ 
of the full harmonic order n < N. In Section 2 we also define the continuous extension 
of the discrete transform on the continuum [0,To]. We then compare CEDCT and 
CEDFT, and show that although both DCT and DFT formally belong to the group 
of exact discrete transforms, surprisingly, only the CEDCT converges with increasing 
A'^ to the continuous function g{t G [0,To]) which originates the grid function {gk}- 

In Section 3, we prove some important properties of the CEDCT. These proper- 
ties closely resemble those of the canonical CFT polynomials where the coefficients 
are found by accurate integrations, such as the principle of locality of CEDCT. This 
feature can ensure, in particular, that numerical computation errors or uncertainties 
in one segment of the interval [0, Tq] would not dramatically affect the results in a 
distant segment, which maybe very important for effective truncation of the discrete 
transform sequence resulting in the loss of the exactness of the discrete transform. 
Another important property proven in Section 3 is that, similar to the CFT polyno- 
mials, the term-by-term derivative series of the continuous extension of an A^-interval 
DCT converges to g'{t) with the increase of A^. Note that this property holds for any 
smooth originating function g{t), in particular when go ^ gp^ (which is not necessarily 
the case for other types of discrete Fourier transforms). 

In Section 4, we extend the formalism of one-dimensional DGT on SU{2), or the 
DCT, for decomposition of multidimensional functions, and bring examples of approx- 
imation of some 2-dimensional discrete functions/images by a continuous extension 
of 2-dimensional CEDCTs. 

2 Basics of Fourier analysis on SU(2) 

The Lie group SU{2) can be realized as a set of all 2x2 unitary complex- valued 
matrices A, with det A = 1. A complex valued class function f on SU{2) is any map 
of SU{2) onto the complex number space C which is invariant under conjugation, 
i.e. / :SU(2) — > C, and f{B-^AB) = f{A) for aU A,B eSU(2). Since the defining 
2-dimensional representation of SU{2) is faithful, we can use it in order to describe 
the discrete elements of SU (2) of interest to us. 

Any unitary matrix can be diagonalizcd by a unitary transformation. Therefore 
every element of SU{2) is conjugate to at least one diagonal matrix in the defining 
2-dimensional representation. All the elements which can be simultaneously diagonal- 
izcd form a maximal torus T of SU{2). Since all maximal tori are S'[/(2)-conjugate, 
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we can write: 



T={x(e)=("7 J^,,^ I 0<^^<l} . (1) 

Furthermore, using B = (i'lj) gSU(2), we have Bx{9)B^^ = x{—6). Therefore 
every element of SU{2) is conjugate to just one element in the subset F C T, where 

Trace functions, otherwise called characters, play an important role in SU{2). For 
any element x{9) of ((T)), we have trx(^^) = 2 cos(27r^). In general, trx(^) is a class 
function because for all B eSU(2)^ tr {B-^xB} = tix. Let R = i?(SU(2)) denote 
a complex algebra generated by the character functions of SU{2). It is well known 
then that R has a linear basis consisting of the characters of all finite dimensional 
irreducible representations of SU{2). 

With each irreducible representation of a semisimple Lie group, in particular SU{2), 
one associates a set of weights (weight system) of the representation [H [3] , which is 
a union of orbits of weights under the action of the Weyl group, W. In physics the 
SU (2)-weights are known as projections of the angular momenta which have integer 
and half-integer eigenvalues. In mathematics one usually prefers to deal with the 
doubles of the angular momenta in order to avoid non-integers. 

The Weyl group of SU{2) is very simple. It is of order 2, consisting of 2 elements 
generated by action of the reflection operator f. It acts on any element m G M of the 
1-dimensional space M of the 'projections of angular momenta' in the straightforward 
way: f(m) = — m. The finite dimensional irreducible representations of SU{2) and of 
its Lie algebra su{2) are well known, but we would point out the following. The 'an- 
gular momentum states', the basis vectors of representation spaces, are eigenvectors 
of the 'diagonal' generator of su{2). Unlike the common normalization of that gener- 
ator in physics, we normalize it so that its eigenvalues are twice the usual projections 
of angular momenta. The set of the eigenvalues designates the weight system of 
the representation /. The weight system of an irreducible representation consists of 
/ + 1 weights, 

n{l) = {mcu I m G {-/, -/ + 2, . . . , / - 2, /} } , 

where / is the highest weight ('twice the angular momentum') of the representation^; 
/ is used to specify the representation. Note that all the elements of have the 
same parity. 

A W^-orbit of a weight m thus consists of one or two elements: 



Wm 



{m, —m} for m ^ 
{0} for m = 



^For sets of weights in case of Lie groups different from SU{2) see [3] 
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The character x{^) is a function of conjugacy classes of the elements of SU{2). 
Every class is represented by one value of 6 within < 6 < ^. The values of the 
character of an irreducible representation / can be written as the sum of values of the 
orbit functions $^(6*), 

x.(o) = t = *,(») + +•■■+{ ; "^t . (3) 

m ^ 

^rn{d) = < , P ^ , and Q < 6 < - . (4) 

^ ^ I 1 for m = ' ~ ~ 2 

Only for the 1- and 2-dimensional representations, / = and 1 respectively, the 
character consists of a single orbit function. Note that $m(6') is symmetric (antisym- 
metric) with respect to the midpoint of its range of 9 for m even (odd). 

The decomposition of irreducible characters Q into the sum of orbit functions is 
given by a triangular matrix. Hence it is invertible. Therefore the orbit functions 
{^m{9)i 1^ = 0,1,2,...} also form a basis in the space of class functions f{6) of 
SU{2). Using (jlD, it implies that 

oo oo ^ 

f{e) = Y,(^m^m{e) =ao + 2^a^ cos(27rme), < < - . (5) 

m=0 m=l 

There are two properties of the expansion © which we want to underline: 

(i) It can be reduced to the familiar case of the standard Fourier decomposition of 
f{9) in the interval 9 G [—1/2, 1/2], if one makes an even extension f{9) = f{~9) for 

te [-1/2,0]. 

(ii) Although f{9) is being expanded into a series of functions which are periodic 
within the range < ^ < 1, the actual range of 9 in ^ makes periodic only the 
cosines with even values of m, i.e m = 2k. Their arguments vary over the range 
{0,2A;7r}, i.e. over an integer multiple of 2tt. 

2.1 Discrete Fourier transform on SU{2) 

The discrete Fourier transform differs from ^ by the fact that the independent 
variable 9 takes only finite number of rational values within its range of variation. 

Fixing a rational value of 9, one fixes an element of finite order (EFO) belonging 
to the SU{2) torus T. Every conjugacy class of EFO in SU{2) is represented by an 
element of T with < ^ < i. In SU{2) one can be explicit, see ^2] or §4 for all 
other compact simple Lie groups. 

Let Tjv denote the set of all elements of T whose adjoint order divides A^, where 
is a positive integer. The adjoint order is the order of the element represented 
by matrices of irreducible representations of SU{2) of odd dimensions (i.e. / even). 
There are exactly (A^ + 1) S'f/(2)-conjugacy classes of such elements. Taking the 
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unique diagonal matrix as representative of each conjugacy class, in representations 
of dimensions 2 and 3, we have the following set of matrices 

T^'' = "IxAr^fe = ^ ^ g- w j ^ k = 0,1, . . . , N 

( / \ 

Tg^= {xN,k= 1 II k = 0,l,...,N 
[ V e-^ 

The trace functions of each of these matrices, of size (/ + 1) x (/ + 1), represent the 
characters Xz(^) of representation /, which can be used for decomposition of the 
class functions on SU{2) (and generally on compact simple Lie groups). 

A more suitable basis for such a decomposition consists of orbit functions jHEI- It 
makes possible {practical) decomposition of class functions on groups of high rank, 
e.g. such as Eg jHl [13 CH] • In the case of SU{2), the orbit functions are also much 
closer to the familiar set of exponentials exp(27rim6') used in the standard Fourier 
analysis. 

Note that the elements xjsi^k are equidistant over the fundamental region F = {5 G 
[0,1/2]} of the Weyl group W. The number of elements of the W^-conjugacy class 
m = k/N is denoted by CN,k and is given by 

_ / 1 if A; = 0,iV 
^^''^ - \ 2 otherwise " 

The following definition of a sesquilinear form (/, g) n in the space R of class func- 
tions / and g on SU (2) is a crucial step for our method: 

N 

(/, 9)n = CN,kfixN,k) g{xN,k) , (7) 

k=Q 

where the overline stands for complex conjugation. It is known jSlElj and it can be 
verified by direct computation, that the set of orbit functions | = 0, . . . , A^} 
is orthogonal on the discrete equidistant A^-interval grid with respect to this form. 
More precisely, we have: 



Jfcm 



for < fc, m < A^. 



For further convenience and for comparison of the results with the conventional 
Fourier series, let us instead of orbit functions $m(6') defined by (jH) consider the 
functions ipmi^) = cos(27rm6') for all m > 0. Then one can easily verify that for 
< A;,m < A^ 

Ctvj cos — - COS — — = — 5km , (8) 

j=0 
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where CN,k is given by (0)- The method proposed below for decomposition of the 
class functions into series of orbit functions is based on the discrete orthogonality 
relation (jH}. 

Let f{9) G -R be a class function that can be decomposed into a finite series of orbit 
functions: 

^ 1 

f{e) = Y,^m^m{e), (Q<e<-). (9) 

m=0 

This can be compared with the general case of infinite series ^ and with the discrete 
Fourier transform (2) in [T7] . 

Our goal now is to find the expansion coefficients {am}- In order to use the or- 
thogonality property (jH)), we form a system of (A^ + 1) linear equations for {flm}, 
restricting 6' in (jH)) to the discrete set of its values {9k = kA6 \ k = 0, . . . , N}, with 

^(^ = m-- 

N N , 

fk = fjOk) = ami'm{Ok) = Qm COS . (10) 

After multiplication of ()10|) by C^^k'^'jk and summing over k, we arrive on the right 
hand side at '^^=Q0,m{4'm,i^j)N- Then, given (jH}, we find 



TV 



J2K'fk forj = 0,...,A^, (11) 



k=0 



where 



j^jk _ Cnj CN,k Tiki 
^N- 2 AT °' AT ■ ^ ' 

Here are the elements of (A^ + l)x(A^+ 1) matrix Dn of the DGT on SU(2). Note 
that it is easily reduced to the transform matrix of the discrete cosine transform of 
the type DCT-1 after renormalization by a factor ^y2CN,j CN,k/N. The matrix Dn 
is independent of the values {fk} of the class function which is being decomposed. 
It is therefore possible to compute in advance, for any predefined values of the 
positive integer A^, and use it repeatedly whenever it is needed. 
Examples of the transform matrices D ^ for the lowest values of A^ are the following: 




--3 



1 1 1 i\ 

1 V2 -V2 -1 
10-20 1 
1 ~V2 V2 -1 



Introducing single-column matrices = {aj} and F/v = {/Ar(ifc)}, and a square 
matrix \I/Ar = {cos(7rjA;/A^) \ j, k = 0, . . . , N}, equations (fTUj) and (fTTj) are written in 
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the matrix form, 



Fn = "^nAn and = D^Fn , where = ■ (14) 

The matrix D^, being the inverse of \E'7v, formally solves the problem of discrete 
Fourier transform on SU{2). 

2.2 Continuous extension of the discrete Fourier transforms 

Let a continuous function g{t) be the origin for the discrete function {g^ = g{tk)} 
defined at the (A^ + 1) points = kTo/N, k = 0,1, . . . N, of the interval [0, To]. The 
DCT of {gk} with the use of the transform matrix |T2|l results in the discrete function 
{aj} in the frequency space. This is an exact ('lossless') discrete transform, since it 
allows unambiguous recovery of all + 1 values of {gk} by applying the inverse DCT 
in the form of (fTUI) . We recall that the standard discrete Fourier transform, i.e. DFT, 
has the same property. 

It seems then natural to ask if it is possible to recover the originating function g{t) 
by a Fourier series not only at the grid points {tk}, but also on the entire continuous 
segment [0,To]. In order to answer this question, let us consider the continuous 
extension of the discrete transform between the grid points, which can be formulated 
as follows. 

Proposition: Let g{t) be a complex valued function of t E [0,To], taking values 
Qk = gitk) on an equidistant point grid {tk = kT^/N \ k = 0,1, . . . , N}. The function 

N 

fN{t e [0, To]) = ^ a, cos(7rjt/To) (15) 

j=0 

with the discrete transform coefficients 

ECN,jCN,k T^jk , s 

-—-—g.COS—, (16) 

A;=0 

represents a continuous Fourier extension of the inverse DCT of {gk}, which is exact 
in the sense that f{tk) = gk all N + 1 points of the grid. 

The proof of this proposition is obvious if one recalls that (fT^ and (fT^ are reduced, 
respectively, to (fTUj) and (fTTj) after substitution t = 29Tq with 9 G [0, |]. 

Similar to DCT, one can also continuously extend to all points of the segment [0, To] 
any other types of discrete transforms, in particular the DFT, resulting in CEDFT. 

Below we will address two important questions. First, how well the CEDCT ap- 
proximates any 'reasonably behaved' (e.g. continuous) function g[t) on the interval 
[0, Tq] outside the points of the grid. The second question is how it compares with 
the CEDFT. At last, we will also briefly address the question of the possible use of 
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CEDCT in practical applications, in particular, for purposes of smooth representation 
of compressed images. 

In order to provide a partial answer to the first question, we consider two examples. 
Example 1. Let us take a Gaussian function 



(7(t)=exp[--(t/(T)2], te[0,l], 



(17) 



with the dispersion a = 1/3, and choose = 3. Thus we have chosen a rather 
coarse grid relative to the dispersion: the width of its intervals 1/A^ is equal to the 
dispersion. The coefficients (|TT|l are readily calculated using from (fT!?j) : 



A, 



ai 



/0.5 + e-°-5 + e-2 + 0.5e-^-5\ 



1 + e 



'0.5 
,-0.5 



e~2 - 



,-2 



,-4.5 



\0.5 + e-°-5 + e-2 + 0.5e-^-7 



/0.415807\ 
0.486695 
0.089748 

\0.007750/ 



(18) 



The corresponding CEDCT function (fT3j) reads: 

fs(t) = ao + ai cos(7rt) + 02 cos(27rt) + 03 cos(37rt). 



(19) 



where the coefficients are given by (fTHj) . 

Now we are in a position to compare (fT7|) with (fT^. At the grid points tk = k/3 the 
functions coincide, fsitk) = g{tk) (which is easily verified). One may then expect large 
deviations of f^it) from g(t) around intermediate points of the grid, i.e. t = 1/6, 1/2, 
and 5/6 : 



73(l/6)\ 
/3(l/2) 
j3(5/6)/ 



^0.88217r 
0.326059 
,0.039191. 



^^(5/6)\ /0.882497\ 
^(1/2) = 0.324652 . 
,^7(1/6)/ VO-043937/ 



A good agreement between /3(t) and g{t) is apparent even for t 
becomes very small. 



5/6 where g{t) 



Example 2. DCT approximation for a more complicated function is found in Fig. 1, 
where we show /Ar(t), = 10 and 14, for the function g(t) composed of two Gaussians, 



git) 
2, A2 



Aie 



_i( tzlz 

+ A2e H -2 . 



1.5, narrow dispersions cti = (J2 



(20) 

0.05, and centered 



with amplitudes Ai 
at ti = 0.42, t2 = 0.56. 

For comparison, in Fig. 1 we show by dashed lines the approximations to g{t) 
provided by the trigonometric CFT polynomials where the transform coefficients are 
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calculated by exact integrations. Recall that for a real function g ^ W the CFT 
polynomials of the harmonic order K are given by the series^ (see e.g. 0): 



K 



K 



PK{t) = J2 



.i2-Kjt/To 



Co + 2 Re Cj 



.i2njt/To 



(21) 



j=-K 



where 




g{t) e-'2-i*/To _ 



(22) 



In order to compare approximations to g{t) that can be provided by series ()15|) and 
{]'21\i with the same order for the highest harmonics, for the series Pxit) in Fig.l we 
put K = N/2. 

Figure 1 illustrates that the DCT is really an exact discrete Fourier transform, i.e. 
that fN{tj) = g{tj) for all < j < + 1. Remarkably, its continuous extension 
/Ar(t G [0,To]) approximates g{t) practically as well as the accurate CFT trigono- 
metric Fourier series Pxit) of the same harmonic order, i.e. K = N/2. In the case 
of Gaussian-type functions, CEDGT series approximates the original function g{t) 
reasonably well even in the case of narrow structures with dispersions as small as 
a ^ To/1.5N. 

2.3 Comparing with standard DFT 

First, let us recall some properties of the standard exploitation of the discrete Fourier 
transform. Further details can be found in many books (e.g. see jSlElEl)- 

The standard DFT is formally derived from an approximate calculation of the 
integral coefficients cj for the trigonometric Fourier series, using a simple rule of 
rectangles for integration of Cj in fl22|) when the function g is given on the A^-interval 
equidistant grid, g{t) {g/, | A; = 0, . . . A^}. This leads to DFT coefficients 



where 1/A^ is the length of the sampling interval At (assuming for simplicity Tq = 1). 

A crucial feature of this definition of DFT consists in the fact that the system of 
equations (f^^ for the first A^ coefficients {uj | j = 0, . . . , A^ — 1} can be inverted with 
respect to {gk \ k = 0, . . . , N — 1}. Such a possibility is based on the observation that 
the matrices with matrix elements 




N-l 



(23) 



k=0 



Mjk = 4^"'^"' and (M-i),fc = e 



^In general, for a complex function g{t) the coefficients C-j ^ cj. 
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are inverse to each other. Thus one gets the inverse DFT in the following form 



N-l 



9k = Y^ Uj 

3=0 




(24) 



Thus, the discrete sets {uj} and {gk} represent a pair of exact (lossless) direct and 
inverse transforms (e.g. see JU]) in the form of Fourier series, which is generally 
treated as the exact solution to the problem of Fourier transform of discrete functions 
given on the equidistant grid. 

We argue below, however, that there are significant reasons to suggest that the real 
(and not only formal) exact solution to the problem of discrete Fourier transform of 
a grid function is provided by the DCT transform pair given by (fTBj) and (fTUI) . It 
retains all the 'good' properties of the DFT: 

(i) it is easy and fast to compute ; 

(ii) it is a lossless discrete transform, with the exact inverse DCT at all (A^ + 1) 
points of the grid (even ii Qq ^ gjq, unlike for DFT); 

However, in addition to these properties, the continuous extension of DCT, /Ar(t) 
in (fT3j) . converges to the originating continuous function g{t) with increasing iV, as 
illustrated in the previous section, and as it is proved in the next Section. Only the 
continuous extension of the DCT, and not the DFT, reveals properties characteristic 
to the canonical continuous Fourier transform series. 

It is worth noting here that the very good convergence properties seem to be a 
common feature for the discrete Lie group transforms, as we also demonstrate on the 
example of SU(3) group in the accompanying paper. The basic mathematics of DGT 
has been formulated [^13j for any dimension n < oo. In fact there are as many different 
variants of the method one could use, as there are different semisimple compact Lie 
groups of rank ra, and then within each variant the choice of the points of the grid 
is also far from unique, except for the lowest cases like SU{2). It then provides an 
opportunity to make a choice of appropriate DGT in situations where the choice of 
symmetry is dictated by the experimental data. 

The absence of the convergence property for the continuous extension of the DFT 
given by (|24|) is not easy to anticipate because CEDFT looks like a Fourier polynomial, 
or a cut-off of an 'ordinary' Fourier expansion: 



Similar to CEDCT, it satisfies the equality hi\]{tk) = Qk on the grid points = k/N 
for^ all k < N — 1 . The fact that their continuous extensions CEDGT and CEDFT 
behave quite differently is illustrated in Example 3. It has rarely been emphasized 
that /iAr(t) does not approximate the initial function g{t) between the grid points, 

^Note that the last k = N grid knot can also be included in DFT provided that g{Q) — g{To). 



N-l 




(25) 



j=0 



11 



and it does not converge at all to any continuous function (except for a trivial case 
of g{t) = const) with increasing A^. It is worth citing in this regard [HI, P- 87: "the 
DFT is a Fourier representation of finite length sequence which is itself a sequence 
rather than a continuous function" . 

Example 3. Let h^it) be the CEDFT of N-interval grid function arising from a 
sampling of the continuous function 

g{t) = Ae-C-^f , (26) 

with A = 2 and parameter a = 0.15. Here we have chosen a large value, 6, for the 
power in the exponent in order to illustrate a case with gradients significantly larger 
than in the case of Gaussian functions. 

Solid hues in Fig. 2 correspond to continuous DFT extensions h^it), and the dashed 
lines show the DOT extension /Ar(t). Although hjyit) passes through all ai t = 
tk (shown as full dots), similar to fN(t), its behaviour in between shows profound 
oscillations due to the presence in (^3)) of high-frequency Fourier components, uj = 
27rj/To (if To 7^ 1), with values of y < j < comparable to A^. These oscillations in 
CEDFT do not decrease with increasing A^, but they quickly disappear in CEDCT. 

This behaviour is explained by the fact that at any large A^ the j-th order harmonic, 
exp(— 2z7rjt), strongly varies and changes its sign in a narrow interval At = 1/A^ for 
high orders j > A^/2. Therefore the rectangular integration rule in (f^^ cannot provide 
any reasonable similarity between the canonical CFT coefficients Cj in and their 
standard DFT 'approximations' Uj in ((221) • Effectively, the 'fine tuning' between the 
coefficients for high-order harmonics intrinsic to the continuous Fourier transform is 
lost. 

2.4 Comments on Fast Fourier Transform 

There is a number of ways to make a significantly more accurate approximation of the 
high-frequency coefficients Cj for a function given on the discrete grid. Despite this 
fact, the standard version of DFT defined by ()2H|1 and for the direct and inverse 
discrete transforms has been widely used since the pioneering paper by Cooley and 
Tukey where the first algorithms for fast calculation of this DFT were devel- 
oped. Different algorithms for such fast Fourier transform (EFT) computations allow 
an increase in the speed of practical calculations of DFTs by one or two orders of mag- 
nitude. Thus, a direct 'head-on' algorithm for calculations of {uj | j = 0, . . . , A^ — 1} 
would require about A^^ multiplications and additions. Meanwhile, for special values 
of A^, EFT algorithms can significantly reduce the required number of elementary 
operations, e.g. down to ~ A^log2 A^ in case of A^ = 2". 

The discussion of various EFT algorithms, extensively developed later on by many 
authors (e.g. [12112011211122]), is outside the scope of this article. Here we note only 
that the EFT methods are fundamentally exploiting the property of the standard 
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DFT coefficients that in their complex-value representation i\2'S\i they can be reduced 
to a power-law series of a single complex element Wn = exp{—i2TT/N) as Uj = 
^^ghW}^. This would not be possible if more accurate integration methods to 
deal with the high-order harmonics were used. It should be noted, however, that the 
DCT does allow for the application of FFT methods since it can be formally reduced to 
2A^-point DFT, and a number of efficient FFT-based algorithms have been developed 
for DCT (e.g. [23 EH ESI EE] ) • Moreover, a number of efficient algorithms competing 
with FFT and specific to DCT have been developed (see [TT] for details). 

A simple modification could significantly improve behaviours of continuous exten- 
sions of discrete Fourier transform series based on the use of standard DFT coefficients 
()23j) ■ but not without penalty. Namely, recall that for a sufficiently smooth function 
g{t) in (j221)j the rectangular integration rule provides good accuracy for approximat- 
ing Cj for low-order harmonics, j < N/2. Since FFT algorithms allow fast calculation 
of {uj}, and as far as g{t) can be generally well approximated with the CFT trigono- 
metric polynomials (PT|) with K < N/2 (as shown below), one can only use the ffist 
half of the standard DFT coefficients Uj for construction of a continuous extension 
of the inverse discrete transform sequence. It will then be similar to the series (j^lj) 
truncated to harmonics of order j < N/2, i.e. it represents a series similar in structure 
to the CEDFT sequence fl24|) . but where the high order harmonics are eliminated: 

K 

s^(t) =uo + 2Re J2 «i e'^"^*^^ ' (27) 
i=i 

where the coefficients Uj are defined by (j^Hj) . Note the difference in the multiplication 
factor 2 in this series at j > 1 as compared with the standard DFT extension /lAr(t) 
of 

In Fig. 3 we show that the function sxit) with K = N/2 (dashed curves) can indeed 
approximate g{t) practically as well as the DCT extension /Ar(t). But in this case the 
penalty is a loss of the 'exactness' property for such modified DFT sequence. That 
is, Sxit) 7^ gk at every t = for A; < iV-|- 1. Therefore (P7|) cannot represent an exact 
solution to the problem of discrete Fourier transform. Meanwhile, the series fN{t) 
both satisfies that condition, and rapidly converges to g{t) with increasing N. 

For comparison, we also show in Fig. 3 (dot-dashed line) the function sxit) cal- 
culated for K = 6 < N/2 = 8. In this case the approximation errors [sK{t) — g{t)] 
are larger than the ones when K = 8, because the order of high frequency harmonics 
becomes important for the approximation of features with a dispersion a < T/2K. 

3 Localization and differentiability of CEDCT 

In this section we prove the properties localization and differentiability of the CEDCT 
which are analogous to the properties of the canonical CFT polynomials ()21|) (e.g. 
see |2]). 
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Derive first a useful formula for an N-interval CEDGT ()15p of a grid function {gk}- 
Using (fTBj). and assuming for simplicity To = 1, i.e. t G [0, 1], (fT3j) is reduced to 



^ N N 

fN{t) = CN,kCN,j cos(7rjt) cos{Tcjtk) Qk , (28) 

where tk = k/N. Using 2 cos a cos P = cos{a + /5) + cos(a — /3), the terms containing 
index j can be re-written as a sum of two geometric series: 

N /N-1 N-1 

Y,Cn,j cos{7i jt) cos(7rjtfc) = (-l)*^cos(7rArt)-l + Re e^"(*+*^)^' + ^ e^" 
j=0 Vj=o j=0 

Summing up the series, one gets a compact expression for CEDGT: 

/iv t = E ' 1 -Qk =J2Mkit)9k, for tG 0,1. 

^-^ 2N cosvrti. — cosTrt 

fc=0 fc=0 

(29) 

The values of /jv (t) for any t ^ tn = n/N are found by applying the Lipshitz rule for 
the ratio of infinitesimals of a smooth function. For (j29p it results in fN{tn) = On for 
all = 0, . . . , A^, as expected. 

Using ()29p. below we prove that the localization principle, known (see j2j) for the 
CFT series ()21|) also holds for the CEDGT. The localization lemma can be formulated 
as follows: 

Localization Lemma. Let the set of N-interval grid functions {gk}, with various 
N , be originated from a smooth function g{t) with a bounded derivative g'{t) on the 
interval [0, 1]. Then at any given t G [0, 1] and for any small A > and e > there 
exists Ai"e,A such that for all N > N^^a the behaviour of the continuous extension of 
the discrete group transform /7v(t) is defined within the accuracy e only by the values 
of g{t) in the /^.-neighborhood oft, i.e.: 
for any e, A > SA^'^^a such that for all N > iV^^A =^ 

Ek' 
AN,k'{t)gk' |< e , with all k' within t — A< — <t + A. 



{k'} 



Proof: From (|T6|) . using (jH}, it follows that the DGT coefficients of any constant 
function gi{t) = const are equal to a^^^ = 6ko x const, i.e. that all coefficients except 
for Oq^^ = const are equal to 0. Thus the convergence properties of the CEDGT of any 
function g(t) are the same as the properties of the function g(t) + const. So taking 
into account that the function is bounded, it is sufficient to prove the Lemma, while 
generally assuming that g(t) is positive on the interval [0,1]. 
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Let us split the sum in into 3 parts corresponding to 

9k , 

k=0 

N 

S+ = ^ AN,k 9k , (30) 

k=Kt+Kc^+l 

Sa = AN,k9k ■ 

k=Kt~KA 

where Kt = [Nt] and Ka = [NA] are the integer parts of the respective products. If 
any of the points t ± A is outside the interval [0,1], then only 2 sub-series are left. In 
all cases, only Sa of these sub-sums is defined by g(t) in the close neighborhood of 
t. The localization Lemma for ()29p then implies that the non- /oca/ sums, 5"+ and S^, 
are reduced with increasing N to absolute values below any small e. 

Indeed, on the basis of (jSHI), for any fixed A each of those non-local sums represents 
a series of bound-value elements with alternating signs, which can be then combined 
into pairs of consequtive elements of order 0(l/iV^) each. Considering for example 
the sum S^, and using the Taylor series decomposition for g^+i = gitk^i) = Qk + 
g'(tk)/N + o{l/N) each pair of elements AN^k9k + AN,k+i9k+i starting with even k > 2 
can be reduced to 



sinvrA^t sinvrt 

iV2 



9'{tk) 9ksimitk ^ o{l/N^ 



cos ntk — cos irt (cos vrt^ — cos irty 



The expression in the square brackets is bounded with some absolute value indepen- 
dent of N as far as g' is bounded and \ tk — t \> A, therefore A^^kgk + ^7v,fc+i5'A;+i ~ 
0(l/iV^). The number of such pairs in S- or S+ is increasing oc A^. Therefore taking 
also into account that in the sums and 5*+ each of the limited number of elements 
left out from such pairing process (e.g. in case of C^^k = 1 for /c = 0, A^) is only of 
order 0(1/A^), we conclude that both series or 5*+ tend to zero with the increase 
of N . This proves the localization Lemma. ■ 

A very important property of the CEDGT series /Ar(t) is the possibility of its term- 
by-term differentiation such that the resulting series converges with increasing A^ to 
the derivative g'{t). Note that while being well-known for the CFT series ()21|) . this 
property is not trivial for finite (N-)element discrete Fourier transforms. Recall that 
although in trigonometric polynomials in the form of (fT3jl or (jTfj) each individual 
term, being oc 1/A^, decreases to at large N, their derivatives over t corresponding 
to a high-order harmonics, say j > N/2, become of order j/N ~ 1, and therefore 
might not necessarily vanish with increasing A^. Thus, a 'fine tuning' of the entire 
discrete-transform based series is needed in order to provide convergence of the series 
produced by its term-by-term differentiation. 
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Theorem: Let g(t) be a smooth function with bounded second derivative on the in- 
terval t G [0, To], which originates the N-interval grid function {g^ \ k = 0,1, . . . , N}. 
Then the function f'j^{t) produced by the term-by-term differentiation of the contin- 
uous extension of the discrete Fourier transform on SU(2), fN{t), converges with 
increasing N to g'{t) at any t G (0,To). 

Proof. As earlier, we put To = 1 for simplicity of the formulae below. Because 
the series fN{t) given by contains a finite number of elements, it is obvious 
that its derivative f'^ can be summed up to the series produced by the term-by-term 
differentiation of ()29|) . 

Consider first the derivative f'p^{t) at any rational to in the open interval (0,1). In 
that case we can choose N such that to = rn/N, and it then makes sense to choose 
all further subdivisions of the interval [0,1] such that to will be always kept as a knot 
of the grid, i.e. to = rrii/Ni for all A^^i > A^, i.e. choosing Ni = aN with some integer 
a > 1. Using the Lipshitz rule, the derivative of at to = m/A^ is reduced to 

TV (kf^m) 

f'M = ^smnto i-^)'~"'CM,kUk{to) + ^^r^^gm, (31) 

z ^ — ' Z sm vrtn 

fc=0 " 

Ukito) = , where t^ = k/N . (32) 

cos vrtfc — cos vrto 

Let us choose some small A, such that both (to ± A) G (0, 1), and then split the 
sum in fl^ into 3 sub-series S'^, S'_ and 5"^ as in (jHUI) . where the number Kt = m for 
t = to = m/N. It is convenient for further analysis to choose as the maximum 
integer which satisfies the condition K^/N < A and is of the same parity as m. 
Then the indices of both the last element in the series S'_ and the first element in S*^, 
m — K/^ — 1 and m + K/^ + 1 respectively, are odd. 

Recall that for any smooth function f/(t) with a bounded second derivative on the 
equidistant grid we have 

2Uk = Uk^i + Uk+i + U'j^ Ar-2 + o(Ar-2) ^ ^33^ 

which follows from the familiar Taylor series decomposition g{t-\-x) = g{t) + g'{t) x + 
0.5 g"{t) x'^+o{x'^) in the X = ±1/A^ vicinity of t. Applying (1221) to f/(t) = 5'(t)/[cos7rt- 
cos vrto] at all t = t^ in dHD) and ^ with k odd (i.e. k = 2j-l \1 < j < (m-i^A) /2 ), 
and given the series S'_ is reduced to 

. , / im~K^)/2 

S'- = {-ir^^^i-U^-K^ + N-^ X: U^^-i\+oiN'')- (34) 

The first term on the right-hand side is exactly one half of the first term in the 
localized sum S^, but with a negative sign. The second term is of order 0{N^^). A 
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more precise estimate of this term can be derived if we note that for large the sum 



{m-K^)/2 

E ^2,-1-/ U"{t)dt = U'{to-A)-U'{0), (35) 
i=i 

is bounded for any given A, as far as g"(t) is bounded. Note that the estimate of 
imphes only that U"{t) is integrable. This suggests that in principle the conditions 
of the validity of the differentiation Theorem can be relaxed, requiring that g"{t) be 
an integrable function on [0, 1], but not necessarily bounded. 
Applying the same approach to S'^, we find 

S'_ + S'^ = !^^^(-l)™+i(f/„_^^ + U^^kJ + Oi , (36) 

with Oi ~ 0{N~^). Note that in the case of an odd N, the residual Oi also includes 
the difference between the last term, k = N, and one half of the k = {N — l)-th term, 
which is of order U'jy sin^nto) / N . Here we take into account that Cn,n-i = '2Cn,n 
from ®. Thus, for any e we can chose such that for all N > the residual 
in (jHUj) is Oi < e/2. Then (jHT|) is reduced with accuracy < e/2 to a summation of 
elements localized around, and symmetric with respect to, the point = m/N: 

Ka-1 

/^(to) = TTSinTTtp {-iy{U,n-j + Um+j) + 

i=i 

+ f.^-J + + 0,(,/2) . (37) 

Z / sm tczq 

Here we take into account that K/^ is chosen to be of the same parity with m. 
Introducing now rj = j/N which is < A for j < K^, we have 

/rr TT \ cos Trio 2 , 

sm 7rto{Um-j + Um+j) = Qm- 1 9m + , 

sm ttzq tt 

where, keeping the largest order terms, the residual O2J is reduced to 

2 , 2 sin^ vrto //\ f 



0'.,^[9,.--C^^^C)[^J . (38) 

Substituting these 2 relations into pTjl . and recalling that g'm = g'(to = m/N), it is 
easily shown that for both even and odd m and Ka one has: 

fM = g'ito) + Oi(e/2) + O^iA') , (39) 

where the residual O2 represents the sum of residuals 02,j, i.e O2 ^2j- Because 
of the sign alteration term (—1)-' in ()37p. this sum does not increase with increasing 
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N beyond the absolute value of O2J at rj = A in ()38p. It follows then that for any 
small e > we can first choose an interval A proportional to ^/e sinvrto (depending 
also on g, g', g") such that | O2 \< e/2. Then we can choose a number N^, such that 
I Oi \< e/2. Hence \ Oi + O2 \< e for all N > N^. This proves the differentiation 
Theorem. ■ 

Note that in the derivation of (j35|) which has allowed principal clipping of both 
end-terms in the localized (jHTj) exactly by one half, the symmetry properties of the 
SU{2) DGT series expressed in the term CN,k of (jHI) have been fully exploited. 

Another notice is that, along with the continuous DGT extension of {gk}, the 
localization Lemma is also valid for its term-by-term derivative f'jyit). This property 
is actually proven by the localized structure of the sum in the construction 

Example 4. In Fig. 4 we show approximations to the function 

^?(t)=e-^* + VK^r (40) 

provided by the continuous extensions of the discrete Fourier transforms in forms 
of the series fll5j) for SU{2) DGT and the series ()27|1 based on the first j < N/2 
half of the standard DFT coefficients uj of ()2Hj) . It is noteworthy to recall that a 
straightforward use of the continuous extension of the standard DFT given by 
for calculations of derivatives does not make any sense as far as such an extension does 
not converge, as shown in Fig. 2. Although the DFT series SN{t) does converge to g{t) 
with increasing N, as shown on the top right-side panel, its derivative series contains 
profound oscillations at any large N. One could suggest that g'{t) can be in principle 
recovered from s'^(t) after the application of some smoothing procedure. Although 
this might be possible, special care should be excercized in order to avoid accumulation 
of systematic errors in circumstances where the amplitude of the oscillations is much 
higher than the mean expected value. Meanwhile, /jv(^) shown by solid lines on the 
two bottom panels of Fig. 4 already provides a rather good approximation to g'(t) at 
relatively low values of A^. 

At the end of this section we would just like to note, albeit without any proof or 
demonstration in this paper, that our numerical calculations show that the second 
derivative of CEDGT also appears to converge to g"{t) if the condition g{0) = giT^) 
is satisfied. 

4 Multidimensional Fourier transform 

Being a transform with separable variables, the multidimensional DGT is easily re- 
duced to the product of one-dimensional DGTs, which is widely used, for example, 
for effective 2D-image processing (e.g. see fUEj. Although the multidimensional 
DGT is well known, in this section we will first formulate it in terms of discrete 
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Fourier transform on the SU{2) x • ■ ■ x SU{2) group. Then we will briefly discuss the 
convergence properties of the continuous extension of 2D discrete cosine transforms. 

The generalization of the transform formulae for decomposition of functions of n 
variables, G{xi, . . . ,Xn), into the Fourier series of orbit functions of [SU(2)]" group 
is straightforward. Let us consider first the case of = 2, i.e. when a function G —>■ 
G{x, y) defined in the region Fn = {0<a;<l, < y < 1} (i.e. assuming normalized 
variables x x/Xq, y u/Yq) is to be decomposed into the series of the orbit 
functions ^rnn{x,y) of the symmetry group SU(2) xSU(2). In this case ^rnn{.x,y) = 
^m{x) ^niy)- Using for convenience again the functions ipmix) = cos^nmx) instead 
of ^rn{x), we can write 

^mn(x, y) = ipm{x)iJn{y) = cos^Tcmx) cos^Tiny) , (41) 

where (x, y) G F. For a uniform rectangular grid {xj, yk} defined in the region F such 
that 

{Xj=j/M, yk = k/N \j=0,l,...,M; k = 0,l,...,N} 
the functions "^mn are orthogonal in the following bilinear form 

M N 

yk)'^pq{Xj,yk) 

j=0 k=0 

which follows directly from (jH)). 

Let G{x, — ^ M be a continuous function originating a 2-dimensional grid function 
Gjk = G{xj,yk) on the grid {xj,yk}. Then the decomposition of G into the Fourier 
series on SU(2)xSU(2) group corresponds to solving the system of equations 

'^jk = X] ^mn "^mniXj, = ^ ^ Amn COS COS , (43) 

m=0 n=0 m=0 n=0 

with {0<j<M,0<A;< A^}, with respect to the coefficients A^n- This 
can be easily achieved using the orthogonality relation ()42|) . if we multiply PHj) by 
CM,jCN,k'^pq{xj, yk) and take the sum over {j, k}. Thus we find the coefficients Amn 
of the discrete Fourier transform (j43|) . 

M N 

= E E ^M^fG^k . (44) 

j=0 k=0 

The matrix is defined as before by (fT^ . 

T^ab CN,aCN,b T^oh 

Dn = """^ 1^ ' N,a,beZ, 
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with the weights CN,a^ CN,b given by 

In this way the coefficients of the 2- dimensional DGT are found for any grid 
function {Gjk} with bounded values at the grid points {xj,yk) G F. Thus we can 
formulate the following proposition: 

Proposition: Let Gjk = G{xj,yk) be values of a bounded function G{x,y) on the 
rectangular grid points 

x,=jXo/M, yu = kY^/N- j e {0, 1, . . . , M} , e {0, 1, . . . , iV} . 

A trigonometric function given by finite Fourier series 

M N 

Fmn{x, y)=y^Z^ ^mn cos — — COS -— , (45) 
-m=0 n=0 " ^ 

with 

. GM,mGM,jGN,nCN,k ^ Trmj TTnk 

Amn = > > GifeCOS COS , (4d) 

continuously extends the discrete inverse Fourier transform of the grid function {Gjk} 
onto the entire rectangular area (x G [0,Xo],?/ G [OjIq]); Q'^^^ satisfies the equality 

Fmn{x,, yk) = Gjk for all j G {0, 1, . . . , M} , A; G {0, 1, . . . , A^} . 

Furthermore, ifG{x, y) is continuous, then Fmn{x, y) converges to G{x, y) for M, N — 
oo. 

Since for any fixed y^ G [0, Fq] or xq G [0,Xo] the series FMN{x,yo) or FMN{xo,y), 
respectively, are reduced to one- dimensional CEDCT/CEDGT series along the x or 
y axes considered in the previous section, it is obvious that the continuous extension 
FMN{x,y) of 2D DGT on the SU(2)xSU(2) group (i.e. the 2-dimensional DGT) 
not only converges with increasing {M,N) to G{x,y), but also that it has properties 
of locality and differentiability similar to the one- dimensional GEDGT on ordinary 
SU{2) group. 

Example 5. The upper panels in Fig. 5 show the contour plots of a function G(x, y) 
defined in the square region F = [0, 1] x [0, 1] and composed of two 2-dimensional 
Gaussian functions, each of type 

e ^ , (47) 

where a\\ > a±, but with directions of the major axes x[ and perpendicular to 
each other. For both Gaussians we have taken the transverse dispersions to be ai^± = 
c"2,± = 0.025, which is exactly 2 times smaller than the grid's cell size Ax = Ay = 1/20 
for the chosen M = N = 20. The contour plots shown on the upper panels in Fig. 5 
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illustrate that even in the case of a grid with cell size this large compared with ctj,, the 
continuous extension of the two-dimensional DGT series reconstructs Gaussian-fast 
smooth structures reasonably well. 

This is also apparent on the bottom panels of Fig. 5 where we show the 3-dimensional 
images for the same analytic function G{x,y) (left panel) and its approximation in 
the form of 2-dimensional continuous DGT extension FMN{x,y) (right panel). Note 
that any waviness that can be seen in the approximated function would disappear 
from the images had we taken N, M >l/a_[_ for the same functions. 
Example 6. The latter case is chosen in Fig. 6 where the upper panels show, in 
terms of brightness distribution, the original grid function produced by 2 Gaussians 
with (Ti,_L = 02,1. — 0.05 = 1/-/V", for = M = 20, and its reconstruction in the 
form of 2-dimensional continuous DGT extension (on the right). The major axes of 
the ellipsoids are inclined at a small angle (20°) to each other. For comparison, on 
the bottom left panel we show the contour plot of the exact (i.e. originating) func- 
tion G{x,y), and the bottom right panel shows its approximation by 2-dimensional 
CEDGT. It is obvious that the reconstructed CEDGT image not only recovers the 
directions of the ellipsoids and their maxima, but it practically coincides with the 
exact image. Note that the dashed contours on both Fig. 5 and Fig. 6 show a level 
slightly below zero, Fmn = —0.001. 

Generalization of the proposition for an n-dimensional DGT of a function G{xi, . . . ,x 
on [SU(2)]" group is straightforward: 

Proposition: LetGj^,,,j^ = G(xijj, . . . ,Xn,j„ be values of a bounded function G{x^^\ . . 
given on the rectangular (Mi, . . . , Mn) -interval grid points 

{x^'^' = jiXi/Mi, . . . , x"'^" = jnXn/Mn , jfc = 0, 1, . . . , Mfe : A; = 0, 1, . . . , Ti} . 

A continuous function given by finite Fourier series 

Ml Mn 

Fm^...mS^\, ...,Xn)= > ... > A^i...m„ COS -— . . . • COS — ; , (48) 

^1 



mi=0 mn=0 



where 



Ml M„ 



■^m\...mn ^ ] ■ ■ ■ ^ ] 



J 1=0 jn=0 



CMi,miCMi,ji • . ■ . • CM„,mnCMn,jn 



X 



Gil in COS — — . . . • COS — — (49) 

satisfies the equality 

FMi...Mn {xji , . . . , XjJ = Gji...j„ , for aU jk G {0, 1, . . . , M^.} and A; e {0, 1, . . . , n} . 
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Furthermore, if G{xi, . . . ,Xn) is continuous, then Fmi...m„{xi, ... ,Xn) converges to 
G{xi, . . . Xn) with Ml, . . . , Mn — > oo. 

The proof of this proposition is readily obtained by the method of induction on the 
SU(2) factors of the group. 

4.1 An example of CEDCT application to real images 

In order to demonstrate the potential of the above suggested approach of continuous 
extension of the inverse multi-dimensional discrete group transforms for purposes of 
natural interpolation of discrete images between the grid points, as well as for the 
possibility of data compression and smooth representation of the compressed images, 
we have chosen a 56 x 140 pixel fragment of the well known image "Lena" . The original 
fragment shown on Fig. 7a is strongly enlarged ('zoomed') in order to make visible 
the granularity of the image at its resolution limits. The grayscale color coding of the 
fragment contains all 256 intensity levels, from g = (black) to g = 255 (white) 

In Fig. 7b we show the continuous extension of the original image. It is constructed 
by subdividing each of the initial intervals Ax and Ay into 3 subintervals. This 
procedure increases the density of the grid points (pixels) by a factor 3x3 = 9. 
Calculations are done dividing first the initial 56x140 pixel fragment into ten 28x28 
pixel subfragments, then calculating the CEDCT for each of these sub-fragments. 
It allows us to demonstrate, on the next two panels, the effects of CEDCT image 
reconstruction at the block edges after some image compression is done. Because the 
continuous extension of the inverse DCT can formally extend the values of the initial 
intensity distribution function to values FMN(x,y) beyond the limits [0,255] used 
for the intensity coding, we have linearly renormalized the values of Fnn —>■ Fnn G 
[0, 255]. The positive impact of the higher resolution achieved by the use of CEDCT 
in Fig. 7b, as compared with the original image in Fig. 7a, is apparent. 

Note that the harmonic order of the cosine functions cos(7r?7,a;/Xo and cos(TTmy /Yq 
used in CEDCT corresponds to the modes < m,n < 50. In Fig. 7c we show the 
continuous extension of the image obtained after application of the simplest "low-pass 
compression" procedure (see e.g. [Hj), putting the DCT coefficients Amn for all 
high-order modes with either m or n exceeding rimax = 19. This procedure generally 
removes the high-frequency 'noise' from the image, and compresses the image by a 
factor (29/20)^ ~ 2. No visual degradation of the CEDCT image is apparent. It 
is noteworthy that although the exactness property of the transform in Fig. 7c is 
lost, the edges of the 10 individual blocks, or the sub-fragments, cannot be visually 
distinguished. The block edges become noticeable only in Fig. 7d, where we have 
applied again the CEDCT approach for visualization of the image compressed now 
by a factor of 10. The compression in Fig. 7d is effectively reached by keeping in 
CEDCT series only 10% of the cosine terms with the large-value coefficients A^n, 
and discarding all terms with small- value Amn- For Fig. 7b it has corresponded to 
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the assumption Amn ^ if | Amn \< 0.05 Amax, where Amax is the maximum absolute 
value of the coefficients {Amn |} {excluding Aqq). Obviously, the image in Fig. 7d still 
remains smooth and quite recognizable. 

We would like to note here that it has not been our aim in this paper to reach 
the goal of the best possible image compression. We believe, however, that through 
this paper our examples demonstrate the high potential of the developed approach of 
continuous extensions of DCT, and of the DGTs generally, as considered in the next 
Paper II for the SU(3) group, for purposes of practical applications, and in particular 
for image processing and compression. 

5 Summary 

We have shown that: 

1. A discrete Fourier transform of a grid function {gk \ k = 0,1, . . . , N} on the orbit 
functions of Lie groups, abbreviated DGT, in the case of SU (2) is reduced to the well 
known discrete cosine transform, namely to DCT-I, which is a known type of exact 
discrete transforms, like the standard DFT sequence. 

2. The principal difference between these 2 types of discrete Fourier transforms 
consists in the fact that DCT is based on the functions cos(/c7rt/To) corresponding to 
k < N trigonometric harmonics of both integer and half-integer orders, n = k/2 < 
N/2, whereas the DFT utilizes trigonometric functions of integer n only, but extending 
to orders n < N. This results in vital differences in the subsequent properties of DFT 
and DCT (or DGT generally). 

3. If the function g(t) originating {g^} is a continuous function of t G [0,To], then 
the continuous extension of the (inverse) DCT sequence results in the function /7v(t) 
which converges to the original g{t) with increasing at all t. This property does 
not hold for the continuous extension of the standard DFT sequence, which shows 
profound oscillations between the points of the grid. Note that potentially this feature 
implies significantly smaller vulnerability of DCT to the truncation/approximation 
errors in the process of filtering as compared with the standard DFT. Therefore it 
could be the reason for the superior general performance of the DCT compared with 
the DFT (see |5]). 

4. Similar to canonical continuous Fourier transform polynomials with coefficients 
calculated by exact integrations, the CEDCT series satisfies the principle of locality. 
This property insures, in particular, that the computation errors connected, e.g., 
with noise or uncertainties in one segment of the data will not significantly affect the 
reconstructed CEDCT image on the distant segment of data. This property of DCT 
may become important especially in the process of lossy data compression when the 
property of exactness of the discrete transform is not necessarily preserved. 

5. Similar to the canonical CFT, the CEDCT series /Ar(t) can be differentiated 
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term by term, so that for the (first) derivative series hmAr^oo f'Nit) ~^ 9'{t) 
t G (0,To) provided that the second derivative of g(t) is a continuous (or just an 
integrable) function on the interval [0,To]. For CEDCT this property is vahd both 
when g{0) = g{To) and g{0) 7^ giTo). It does not necessarily hold for other types of 
discrete Fourier transforms which might themselves be converging, like Sk{t) in ()27|1 . 
but which produce nontheless a non- converging derivative series, as demonstrated in 
Fig. 4. 

The derivative series satisfies the localization principle along with /Ar(t). 

6. In the case of an ra-dimensional function defined on the knots of a rectangular 
ra-dimensional grid, the DGT Fourier decomposition can be performed using the orbit 
functions of [S'f/(2)]" group. Such Fourier series are effectively reduced to the n-fold 
convolution of one-dimensional DGT on SU(2) alongside n independent (rectangu- 
lar) axes, and therefore they posess nice properties of convergence, localization, and 
differentiability of their continuous DGT extensions similar to the one- dimensional 
CEDCT. 
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Figure 1: Approximation of an analytic function (j^U]) . shown by the dotted hues, with 
the CEDCT function f^it) (sohd hues) for the discrete interval numbers = 10 
and = 14. The big dots show the values of the grid function {gk = g(tk)} for 
{k = 0, . . . , A^}. For comparison we also show by dashed lines the approximation of 
g{t) by the exact CFT polynomials of (PT|) with K = N/2. It gives polynomials of 
the same harmonic maximum order as the ones in the DCT. The parameters in (I^U]) 
are: Al = 2, A2 = 1.5, = ^2 = 0.05, = 0.42, = 0.56. 
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Figure 2: Behaviours of the Fourier series hiy{t) and /Ar(t), given by equations 
and (|T5|) . Sohd hues represent the CEDFT, and dashed hues show the CEDCT. 
Big dots show the values of the grid function {gk \ k = 0, . . . , N} originated from an 
analytic function g{t) which is given by equation ()2fj|l . and is shown by the dot-dashed 
line. 
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Figure 3: Approximation of an analytic function g{t) composed of two Gaussians, 
(1201), with the CEDCT function /7v(t) of (US)) with = 16 {soM Imes) and the 
Fourier series sxit) of ^Tl\ that uses only the first K ~ N/2 coefficients {uj \ j < K} 
of the standard DFT The dashed line is for K = N/2 = 8, and the dot- 

dashed line corresponds to K = 6. The full dots in the upper panel correspond 
to {gk I k = 0,...,A^}. The lower panel shows the corresponding errors of the 
approximation to g{t) by these 3 types of discrete Fourier transforms. The dispersions 
in (PUI) are assumed to be ai = 0.07 and cr2 = 0.2 . 
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Figure 4: The CEDCT series f^ify {solid lines, top panels) and the derivatives fjsfit) 
[solid lines, bottom panels) for an analytic function g{t) of equation ()40p and its 
derivative g'{t) {heavy dot- dashed lines) in the case of = 14 and = 140. For 
comparison, by dashed hues we show the truncated DFT series sxit) of equation (P7jl 
and its derivative s'j^{t) with K = N/2. Note that although sxit) converges with 
^ oo at all t G (0, 1), its derivative s'j^ apparently does not. 
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Figure 5: Contour plots (upper panels) and visual '3D' images (lower panels) of 
an analytic function G{x,y) (right side) and its approximation by 2-dimensional 
SUi2)xSU{2) CEDGT/DCT series FMN{x,y) (left side) with M = = 20. G{x,y) 
is composed of a sum of 2 two-dimensional Gaussian ellipsoids, each in the form of 
(jUj), and with dispersions a^^i = a±^2 = 0.025 = 1/2N. The dashed lines on the left 
bottom panel show the contour level FMN{x,y) = —0.001 . 



31 



Figure 6: Upper panels: The images, in terms of brightness distributions, cor- 
responding to the grid function {Gj^k} given on M = = 20 square grid (on the 
left), and its reconstruction by 2-dimensional CEDCT series Fm.n{,x, y) (on the right). 
Lower panels: Contour plots of the analytic function G{x^ y) originating {Gj,fe} (on 
the left), and of its approximation by FM,N{x.y) (on the right). The original func- 
tion G{x,y) is composed of two-dimensional Gaussian distributions with transverse 
dispersions for both a± — 1/N — 0.05, an angle between long axes of the ellipsoids 
equal to 20°, and a small separation between their peak positions. 
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Figure 7: (a) Zoomed 56 x 140-pixel fragment of the original image "Lena"; (b) 
CEDCT image calculated with the 3x3 higher resolution, i.e. corresponding to 168 x 
420 pixels; For our calculations the image was first subdivided into 10 square blocks of 
28 X 28 size, (c) The picture in the same 168 x 420 pixel representation compressed 
by a factor of 2. The compression is reached by retaining in only the first 
< m,n < 19 DCT coefficients Amn, and discarding all higher-order terms, (d) The 
picture compressed by a factor of 10; the compression is achieved by discarding all 
cosine terms with | Amn |< 0.05 A^ax (see text). 
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